\(\int \frac {a+b x^2+c x^4}{x^6 \sqrt {d-e x} \sqrt {d+e x}} \, dx\) [143]

Optimal result

Integrand size = 35, antiderivative size = 160 \[ \int \frac {a+b x^2+c x^4}{x^6 \sqrt {d-e x} \sqrt {d+e x}} \, dx=-\frac {a \left (d^2-e^2 x^2\right )}{5 d^2 x^5 \sqrt {d-e x} \sqrt {d+e x}}-\frac {\left (5 b d^2+4 a e^2\right ) \left (d^2-e^2 x^2\right )}{15 d^4 x^3 \sqrt {d-e x} \sqrt {d+e x}}-\frac {\left (15 c d^4+10 b d^2 e^2+8 a e^4\right ) \left (d^2-e^2 x^2\right )}{15 d^6 x \sqrt {d-e x} \sqrt {d+e x}} \]

[Out]

-1/5*a*(-e^2*x^2+d^2)/d^2/x^5/(-e*x+d)^(1/2)/(e*x+d)^(1/2)-1/15*(4*a*e^2+5*b*d^2)*(-e^2*x^2+d^2)/d^4/x^3/(-e*x
+d)^(1/2)/(e*x+d)^(1/2)-1/15*(8*a*e^4+10*b*d^2*e^2+15*c*d^4)*(-e^2*x^2+d^2)/d^6/x/(-e*x+d)^(1/2)/(e*x+d)^(1/2)

Rubi [A] (verified)

Time = 0.10 (sec) , antiderivative size = 160, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.114, Rules used = {534, 1279, 464, 270} \[ \int \frac {a+b x^2+c x^4}{x^6 \sqrt {d-e x} \sqrt {d+e x}} \, dx=-\frac {\left (d^2-e^2 x^2\right ) \left (8 a e^4+10 b d^2 e^2+15 c d^4\right )}{15 d^6 x \sqrt {d-e x} \sqrt {d+e x}}-\frac {\left (d^2-e^2 x^2\right ) \left (4 a e^2+5 b d^2\right )}{15 d^4 x^3 \sqrt {d-e x} \sqrt {d+e x}}-\frac {a \left (d^2-e^2 x^2\right )}{5 d^2 x^5 \sqrt {d-e x} \sqrt {d+e x}} \]

[In]

Int[(a + b*x^2 + c*x^4)/(x^6*Sqrt[d - e*x]*Sqrt[d + e*x]),x]

[Out]

-1/5*(a*(d^2 - e^2*x^2))/(d^2*x^5*Sqrt[d - e*x]*Sqrt[d + e*x]) - ((5*b*d^2 + 4*a*e^2)*(d^2 - e^2*x^2))/(15*d^4
*x^3*Sqrt[d - e*x]*Sqrt[d + e*x]) - ((15*c*d^4 + 10*b*d^2*e^2 + 8*a*e^4)*(d^2 - e^2*x^2))/(15*d^6*x*Sqrt[d - e
*x]*Sqrt[d + e*x])

Rule 270

Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[(c*x)^(m + 1)*((a + b*x^n)^(p + 1)/(a*
c*(m + 1))), x] /; FreeQ[{a, b, c, m, n, p}, x] && EqQ[(m + 1)/n + p + 1, 0] && NeQ[m, -1]

Rule 464

Int[((e_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_)), x_Symbol] :> Simp[c*(e*x)^(m +
1)*((a + b*x^n)^(p + 1)/(a*e*(m + 1))), x] + Dist[(a*d*(m + 1) - b*c*(m + n*(p + 1) + 1))/(a*e^n*(m + 1)), In
t[(e*x)^(m + n)*(a + b*x^n)^p, x], x] /; FreeQ[{a, b, c, d, e, p}, x] && NeQ[b*c - a*d, 0] && (IntegerQ[n] ||
GtQ[e, 0]) && ((GtQ[n, 0] && LtQ[m, -1]) || (LtQ[n, 0] && GtQ[m + n, -1])) &&  !ILtQ[p, -1]

Rule 534

Int[(u_.)*((c_) + (d_.)*(x_)^(n_.) + (e_.)*(x_)^(n2_.))^(q_.)*((a1_) + (b1_.)*(x_)^(non2_.))^(p_.)*((a2_) + (b
2_.)*(x_)^(non2_.))^(p_.), x_Symbol] :> Dist[(a1 + b1*x^(n/2))^FracPart[p]*((a2 + b2*x^(n/2))^FracPart[p]/(a1*
a2 + b1*b2*x^n)^FracPart[p]), Int[u*(a1*a2 + b1*b2*x^n)^p*(c + d*x^n + e*x^(2*n))^q, x], x] /; FreeQ[{a1, b1,
a2, b2, c, d, e, n, p, q}, x] && EqQ[non2, n/2] && EqQ[n2, 2*n] && EqQ[a2*b1 + a1*b2, 0]

Rule 1279

Int[((f_.)*(x_))^(m_)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Wit
h[{Qx = PolynomialQuotient[(a + b*x^2 + c*x^4)^p, f*x, x], R = PolynomialRemainder[(a + b*x^2 + c*x^4)^p, f*x,
x]}, Simp[R*(f*x)^(m + 1)*((d + e*x^2)^(q + 1)/(d*f*(m + 1))), x] + Dist[1/(d*f^2*(m + 1)), Int[(f*x)^(m + 2)
*(d + e*x^2)^q*ExpandToSum[d*f*(m + 1)*(Qx/x) - e*R*(m + 2*q + 3), x], x], x]] /; FreeQ[{a, b, c, d, e, f, q},
x] && NeQ[b^2 - 4*a*c, 0] && IGtQ[p, 0] && LtQ[m, -1]

Rubi steps \begin{align*} \text {integral}& = \frac {\sqrt {d^2-e^2 x^2} \int \frac {a+b x^2+c x^4}{x^6 \sqrt {d^2-e^2 x^2}} \, dx}{\sqrt {d-e x} \sqrt {d+e x}} \\ & = -\frac {a \left (d^2-e^2 x^2\right )}{5 d^2 x^5 \sqrt {d-e x} \sqrt {d+e x}}-\frac {\sqrt {d^2-e^2 x^2} \int \frac {-5 b d^2-4 a e^2-5 c d^2 x^2}{x^4 \sqrt {d^2-e^2 x^2}} \, dx}{5 d^2 \sqrt {d-e x} \sqrt {d+e x}} \\ & = -\frac {a \left (d^2-e^2 x^2\right )}{5 d^2 x^5 \sqrt {d-e x} \sqrt {d+e x}}-\frac {\left (5 b d^2+4 a e^2\right ) \left (d^2-e^2 x^2\right )}{15 d^4 x^3 \sqrt {d-e x} \sqrt {d+e x}}+\frac {\left (\left (15 c d^4-2 e^2 \left (-5 b d^2-4 a e^2\right )\right ) \sqrt {d^2-e^2 x^2}\right ) \int \frac {1}{x^2 \sqrt {d^2-e^2 x^2}} \, dx}{15 d^4 \sqrt {d-e x} \sqrt {d+e x}} \\ & = -\frac {a \left (d^2-e^2 x^2\right )}{5 d^2 x^5 \sqrt {d-e x} \sqrt {d+e x}}-\frac {\left (5 b d^2+4 a e^2\right ) \left (d^2-e^2 x^2\right )}{15 d^4 x^3 \sqrt {d-e x} \sqrt {d+e x}}-\frac {\left (15 c d^4+10 b d^2 e^2+8 a e^4\right ) \left (d^2-e^2 x^2\right )}{15 d^6 x \sqrt {d-e x} \sqrt {d+e x}} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.17 (sec) , antiderivative size = 87, normalized size of antiderivative = 0.54 \[ \int \frac {a+b x^2+c x^4}{x^6 \sqrt {d-e x} \sqrt {d+e x}} \, dx=-\frac {\sqrt {d-e x} \sqrt {d+e x} \left (15 c d^4 x^4+5 b d^2 x^2 \left (d^2+2 e^2 x^2\right )+a \left (3 d^4+4 d^2 e^2 x^2+8 e^4 x^4\right )\right )}{15 d^6 x^5} \]

[In]

Integrate[(a + b*x^2 + c*x^4)/(x^6*Sqrt[d - e*x]*Sqrt[d + e*x]),x]

[Out]

-1/15*(Sqrt[d - e*x]*Sqrt[d + e*x]*(15*c*d^4*x^4 + 5*b*d^2*x^2*(d^2 + 2*e^2*x^2) + a*(3*d^4 + 4*d^2*e^2*x^2 +
8*e^4*x^4)))/(d^6*x^5)

Maple [A] (verified)

Time = 0.45 (sec) , antiderivative size = 82, normalized size of antiderivative = 0.51

method	result	size
gosper	\(-\frac {\sqrt {e x +d}\, \sqrt {-e x +d}\, \left (8 a \,e^{4} x^{4}+10 b \,d^{2} e^{2} x^{4}+15 c \,d^{4} x^{4}+4 a \,d^{2} e^{2} x^{2}+5 b \,d^{4} x^{2}+3 a \,d^{4}\right )}{15 x^{5} d^{6}}\)	\(82\)
risch	\(-\frac {\sqrt {e x +d}\, \sqrt {-e x +d}\, \left (8 a \,e^{4} x^{4}+10 b \,d^{2} e^{2} x^{4}+15 c \,d^{4} x^{4}+4 a \,d^{2} e^{2} x^{2}+5 b \,d^{4} x^{2}+3 a \,d^{4}\right )}{15 x^{5} d^{6}}\)	\(82\)
default	\(-\frac {\sqrt {-e x +d}\, \sqrt {e x +d}\, \operatorname {csgn}\left (e \right )^{2} \left (8 a \,e^{4} x^{4}+10 b \,d^{2} e^{2} x^{4}+15 c \,d^{4} x^{4}+4 a \,d^{2} e^{2} x^{2}+5 b \,d^{4} x^{2}+3 a \,d^{4}\right )}{15 d^{6} x^{5}}\)	\(86\)

[In]

int((c*x^4+b*x^2+a)/x^6/(-e*x+d)^(1/2)/(e*x+d)^(1/2),x,method=_RETURNVERBOSE)

[Out]

-1/15*(e*x+d)^(1/2)*(-e*x+d)^(1/2)*(8*a*e^4*x^4+10*b*d^2*e^2*x^4+15*c*d^4*x^4+4*a*d^2*e^2*x^2+5*b*d^4*x^2+3*a*
d^4)/x^5/d^6

Fricas [A] (verification not implemented)

none

Time = 0.27 (sec) , antiderivative size = 76, normalized size of antiderivative = 0.48 \[ \int \frac {a+b x^2+c x^4}{x^6 \sqrt {d-e x} \sqrt {d+e x}} \, dx=-\frac {{\left (3 \, a d^{4} + {\left (15 \, c d^{4} + 10 \, b d^{2} e^{2} + 8 \, a e^{4}\right )} x^{4} + {\left (5 \, b d^{4} + 4 \, a d^{2} e^{2}\right )} x^{2}\right )} \sqrt {e x + d} \sqrt {-e x + d}}{15 \, d^{6} x^{5}} \]

[In]

integrate((c*x^4+b*x^2+a)/x^6/(-e*x+d)^(1/2)/(e*x+d)^(1/2),x, algorithm="fricas")

[Out]

-1/15*(3*a*d^4 + (15*c*d^4 + 10*b*d^2*e^2 + 8*a*e^4)*x^4 + (5*b*d^4 + 4*a*d^2*e^2)*x^2)*sqrt(e*x + d)*sqrt(-e*
x + d)/(d^6*x^5)

Sympy [F(-1)]

Timed out. \[ \int \frac {a+b x^2+c x^4}{x^6 \sqrt {d-e x} \sqrt {d+e x}} \, dx=\text {Timed out} \]

[In]

integrate((c*x**4+b*x**2+a)/x**6/(-e*x+d)**(1/2)/(e*x+d)**(1/2),x)

[Out]

Timed out

Maxima [A] (verification not implemented)

none

Time = 0.28 (sec) , antiderivative size = 148, normalized size of antiderivative = 0.92 \[ \int \frac {a+b x^2+c x^4}{x^6 \sqrt {d-e x} \sqrt {d+e x}} \, dx=-\frac {\sqrt {-e^{2} x^{2} + d^{2}} c}{d^{2} x} - \frac {2 \, \sqrt {-e^{2} x^{2} + d^{2}} b e^{2}}{3 \, d^{4} x} - \frac {8 \, \sqrt {-e^{2} x^{2} + d^{2}} a e^{4}}{15 \, d^{6} x} - \frac {\sqrt {-e^{2} x^{2} + d^{2}} b}{3 \, d^{2} x^{3}} - \frac {4 \, \sqrt {-e^{2} x^{2} + d^{2}} a e^{2}}{15 \, d^{4} x^{3}} - \frac {\sqrt {-e^{2} x^{2} + d^{2}} a}{5 \, d^{2} x^{5}} \]

[In]

integrate((c*x^4+b*x^2+a)/x^6/(-e*x+d)^(1/2)/(e*x+d)^(1/2),x, algorithm="maxima")

[Out]

-sqrt(-e^2*x^2 + d^2)*c/(d^2*x) - 2/3*sqrt(-e^2*x^2 + d^2)*b*e^2/(d^4*x) - 8/15*sqrt(-e^2*x^2 + d^2)*a*e^4/(d^
6*x) - 1/3*sqrt(-e^2*x^2 + d^2)*b/(d^2*x^3) - 4/15*sqrt(-e^2*x^2 + d^2)*a*e^2/(d^4*x^3) - 1/5*sqrt(-e^2*x^2 +
d^2)*a/(d^2*x^5)

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1055 vs. \(2 (145) = 290\).

Time = 0.53 (sec) , antiderivative size = 1055, normalized size of antiderivative = 6.59 \[ \int \frac {a+b x^2+c x^4}{x^6 \sqrt {d-e x} \sqrt {d+e x}} \, dx=-\frac {4 \, {\left (15 \, c d^{4} e^{2} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{9} + 15 \, b d^{2} e^{4} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{9} + 15 \, a e^{6} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{9} - 240 \, c d^{4} e^{2} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{7} - 160 \, b d^{2} e^{4} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{7} - 80 \, a e^{6} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{7} + 1440 \, c d^{4} e^{2} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{5} + 800 \, b d^{2} e^{4} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{5} + 928 \, a e^{6} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{5} - 3840 \, c d^{4} e^{2} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{3} - 2560 \, b d^{2} e^{4} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{3} - 1280 \, a e^{6} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{3} + 3840 \, c d^{4} e^{2} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )} + 3840 \, b d^{2} e^{4} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )} + 3840 \, a e^{6} {\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}\right )}}{15 \, {\left ({\left (\frac {\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}{\sqrt {e x + d}} - \frac {\sqrt {e x + d}}{\sqrt {2} \sqrt {d} - \sqrt {-e x + d}}\right )}^{2} - 4\right )}^{5} d^{6} e} \]

[In]

integrate((c*x^4+b*x^2+a)/x^6/(-e*x+d)^(1/2)/(e*x+d)^(1/2),x, algorithm="giac")

[Out]

-4/15*(15*c*d^4*e^2*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(
-e*x + d)))^9 + 15*b*d^2*e^4*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d
) - sqrt(-e*x + d)))^9 + 15*a*e^6*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*s
qrt(d) - sqrt(-e*x + d)))^9 - 240*c*d^4*e^2*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/
(sqrt(2)*sqrt(d) - sqrt(-e*x + d)))^7 - 160*b*d^2*e^4*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt
(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(-e*x + d)))^7 - 80*a*e^6*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) -
sqrt(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(-e*x + d)))^7 + 1440*c*d^4*e^2*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt
(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(-e*x + d)))^5 + 800*b*d^2*e^4*((sqrt(2)*sqrt(d) - sqrt(-e*x
+ d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(-e*x + d)))^5 + 928*a*e^6*((sqrt(2)*sqrt(d) - sqrt
(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(-e*x + d)))^5 - 3840*c*d^4*e^2*((sqrt(2)*sqr
t(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(-e*x + d)))^3 - 2560*b*d^2*e^4*((
sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(-e*x + d)))^3 - 1280*a
*e^6*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(-e*x + d)))^3 +
3840*c*d^4*e^2*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) - sqrt(-e*x
+ d))) + 3840*b*d^2*e^4*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) -
sqrt(-e*x + d))) + 3840*a*e^6*((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(
d) - sqrt(-e*x + d))))/((((sqrt(2)*sqrt(d) - sqrt(-e*x + d))/sqrt(e*x + d) - sqrt(e*x + d)/(sqrt(2)*sqrt(d) -
sqrt(-e*x + d)))^2 - 4)^5*d^6*e)

Mupad [B] (verification not implemented)

Time = 8.49 (sec) , antiderivative size = 146, normalized size of antiderivative = 0.91 \[ \int \frac {a+b x^2+c x^4}{x^6 \sqrt {d-e x} \sqrt {d+e x}} \, dx=-\frac {\sqrt {d-e\,x}\,\left (\frac {a}{5\,d}+\frac {x^4\,\left (15\,c\,d^5+10\,b\,d^3\,e^2+8\,a\,d\,e^4\right )}{15\,d^6}+\frac {x^5\,\left (15\,c\,d^4\,e+10\,b\,d^2\,e^3+8\,a\,e^5\right )}{15\,d^6}+\frac {x^2\,\left (5\,b\,d^5+4\,a\,d^3\,e^2\right )}{15\,d^6}+\frac {x^3\,\left (5\,b\,d^4\,e+4\,a\,d^2\,e^3\right )}{15\,d^6}+\frac {a\,e\,x}{5\,d^2}\right )}{x^5\,\sqrt {d+e\,x}} \]

[In]

int((a + b*x^2 + c*x^4)/(x^6*(d + e*x)^(1/2)*(d - e*x)^(1/2)),x)

[Out]

-((d - e*x)^(1/2)*(a/(5*d) + (x^4*(15*c*d^5 + 10*b*d^3*e^2 + 8*a*d*e^4))/(15*d^6) + (x^5*(8*a*e^5 + 10*b*d^2*e
^3 + 15*c*d^4*e))/(15*d^6) + (x^2*(5*b*d^5 + 4*a*d^3*e^2))/(15*d^6) + (x^3*(4*a*d^2*e^3 + 5*b*d^4*e))/(15*d^6)
+ (a*e*x)/(5*d^2)))/(x^5*(d + e*x)^(1/2))